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1. Introduction 

In this report we summarize new results on lepton-nucleon deep inelastic 
scattering (DIS) from the fixed target experiments NuTeV and Hermes, and 
from the ep collider HERA in the range of four- momentum transfer squared, 
Q 2 , between 0.35 and 30000 GeV 2 . Recent next-to-leading-order (NLO) 
QCD fits and uncertainties of the parton distribution functions (PDFs) are 
discussed. In the theory part of this summary, developments in lattice QCD, 
saturation type effects and colour glass condensates, A^-factorization, and 
collinear factorization are surveyed. 
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2. DIS results on the proton structure functions 

2.1. The first NuTeV results on Fi 

The i/(F)-nucleon cross sections from the NuTeV neutrino experiment 
were presented by R. Bernstein [1]. In Fig.l (left) the results are compared 
with the former CCFR data [2] as function of the inelasticity y. The mea- 
surements are in a good agreement apart from x = 0.45 where the v results 
of CCFR are systematically lower. In contrast to CCFR the NuTeV experi- 
ment uses very clean v and V beams provided by a Sign Selected Quadrupolc 
Train where the charge of the parent n, K of the neutrinos can be selected. 
The admixtures of the wrong neutrino type is 3 • 10~ 4 for v and 4 • 10~ 3 for 
V beams. The energy scale uncertainties for muons and hadrons are also 
improved compared to CCFR with 0.8% for muon (goal 0.3%) and 0.4% for 
hadrons. For CCFR both uncertainties were 1%. 

The sum of v and V cross sections depends on F2, R, the ratio o-l/o~t of 
longitudinal to transverse cross sections, and AXF3 = xF^ — xF^ = 4x(s — c) 
which is sensitive to heavy quark densities. All three functions cannot be 
derived from the data simultaneously because of strong correlations among 
corresponding parameters. The first NuTeV results on F2(x,Q 2 ), shown 
in Fig.l (right), were obtained using the world knowledge on R and AxF% 
deduced from the y dependence of the cross sections. 




Fig. 1. v{v)-nucleon cross sections (left) and F 2 (x, Q 2 ) (right) from NuTeV. 
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2.2. New HERMES results on nuclear effects in DIS 

DIS cross section ratios for positrons of 27.5 GeV on helium-3, nitrogen 
and krypton with respect to deuterium were measured by the HERMES 
collaboration (presented by A. Bruell [3]). The helium-3 and nitrogen data 
were already published [4]. Recently, those data were found to suffer from 
an A-dependent tracking efficiency of the HERMES spectrometer, which 
was not recognised in the previous analysis. The resulting correction of the 
cross section ratios is significant at low values of x and Q 2 and substantially 
changes the interpretation of those data. The data corrected for this effect 
are shown in Fig. 2. They are in agreement with previous measurements of 
NMC and SLAC. Values for the ratio Ra/Rd have been derived from the 
the y dependence of the data and are found to be consistent with unity. 
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Fig. 2. Ratio of isoscalar DIS Born cross sections for several nuclei with respect to 
deuterium (left), ratio Ra/Rd (right). 



2.3. HERA results at low and medium Q 2 

T. Lastovicka [5] presented new HI data on ^(cc, Q 2 ) at very low x and 
0.35 < Q 2 < 3.5 GeV 2 in the transition region from the non-perturbative 
QCD to the DIS domain, see Fig. 3. The data were taken in 2000 in a 
special run with the interaction vertex shifted by 70 cm in the proton beam 
direction, thereby accessing lower Q 2 than at the nominal vertex position. 
The luminosity was increased by about a factor of four as compared to the 
initial shifted vertex run in 1995 which lead to the first HI [6] and ZEUS [7] 
data on the proton structure function in the low Q 2 domain. 
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Fig. 3. F2/Q 2 ~ &tot{l*p) as function of Q 2 for different W , the invariant mass of 
the hadronic final state. Grey squares at 0.35 < Q 2 < 3.5 GeV 2 correspond to the 
new HI results. 



The steep rise of the proton structure function F 2 towards small x was 
first observed in 1993 in the HERA data [8]. In perturbative QCD this rise 
corresponds to an increase of the gluon density and is expected to slow down 
at highest energies (small x) due to gluon-gluon interactions. Meanwhile the 
precision of the F 2 data is much improved and the rise is studied in great 
detail. J. Gayler [9] presented the local derivative A = — (d In F 2 /d In x)q2 
based on the new HI F 2 data [5] and published precision HI data [10]. The x 
and Q 2 dependence of A is shown in Fig. 4 (left). The derivative is constant 
for fixed Q 2 in the range x < 0.01 consistent with the QCD fit. Therefore 
the data were fitted assuming the power behaviour F 2 = c(Q 2 )x- x( - Q >. The 
results for the A and c values are presented in Fig. 4 (right). At Q 2 < 2 GeV 2 
the HI data were combined with data of NMC [11] and ZEUS [12]. We can 
state that no damping effects of the rise of F2 are visible yet at present 
energies for Q 2 > 0.85 GeV 2 . For Q 2 > 3.5 GeV 2 and x < 0.01, F 2 can 
be well described by the very simple parameterisation F 2 = cx~ x ^ ' with 
c 0.18 and X(Q 2 ) = adn(Q 2 /A 2 ). At very low Q 2 A is approaching 0.08 
which corresponds to the energy dependence of soft hadronic interactions 
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Fig. 4. Local derivative A = — (d\n F2/dlnx)Q2 (left) and fitted values of X(Q 2 ) 
and c(Q 2 ) (right). 
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24. High Q 2 HERA data 

New high Q 2 HERA data were presented by M. Ellerbrock [14], M. Moritz [15] 
and S. Grijpink [16]. Both ZEUS and HI have results from ~ \6pb~ 1 of e~p 
data taken in the years 1998-1999 and ~ 60p6 _1 of e + p data taken in the 
years 1999-2000, both at y^s = 318GeV. These e + p data can be combined 
with the previously published data at y/s = 300GeV to give a total sample 
of ~ 90P&- 1 . 

The ZEUS and HI data for neutral (NC) and charged current (CC) e^p 
scattering are compared in Fig. 5. There is excellent agreement between 
the experiments and with the Standard Model predictions for electroweak 
unification at high Q 2 . 

Neglecting the small contribution of Fl, the differential cross-section for 
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NC e p scattering is given by 

where F2 , xF 3 are expressed in terms of parton distribution functions (PDFs) 
as 

F 2 (x, Q 2 ) = ZiMQ^ixqiix, Q 2 ) + xqi(x, Q 2 )) (2) 

and 

xF 3 (x, Q 2 ) = ZiBiiQ^ixqiix, Q 2 ) - xq, t {x, Q 2 )) (3) 

in leading order perturbative QCD. For unpolarised lepton beams the co- 
efficients A,B are given in terms of electroweak couplings [18]. The parity 
violating structure function xF% is only significant at high Q 2 . Fig 5 also 
shows the difference in the e + and e~ NC cross-sections due to this xF$ 
term as a function of Q 2 . This has been used to extract xF 3 and a new 
measurement from ZEUS is shown in Fig 6. With the greater luminosity of 
HERA-II a precison measurement of this valence structure function will be 
possible across all x. Currently the only such precision measurement is the 
XF3 measurement from CCFR u, v scattering on an Fe target. 

High Q 2 CC data can be used to gain information on the high x valence 
PDFs, with flavour separation between u v and d v . CC scattering involves 
only the quark flavours which are appropriate to the charge of the current, 
so that the differential cross-section for CC e^p scattering with unpolarized 
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Fig. 6. The structure function xF 3 extracted from ZEUS high Q 2 NC e ± p data. 



beams is given by 



d 2 a(e-p) _ G\ 



dxdQ 2 2irx (Q 2 + M^) 2 



xU(x,Q 2 ) + (l-y) 2 xD(x,Q 2 ) (4) 



and 



d 2 a(e+p) _ G\ 



dxdQ 2 2ttx (Q 2 + M^) 2 



xU(x,Q 2 ) + (1 - y) 2 xD(x,Q 2 ) 



(5) 



where U stands for [/-type quarks with charge +2/3 and D for L>-type 
quarks with charge —1/3. Clearly at high x the e~p cross-section is dom- 
inated by the u v PDF and the e + p cros-section by the d v PDF. The new 
high Q 2 CC data are shown in Fig 7. Their contribution to the precision 
extraction of PDFs will be discussed in Sec. 3.2. 

The strong dependence of the CC cross-sections on the W propagator 
can be used to make an extraction of Myy in a space-like process. The 
results from ZEUS and HI for e~ and e + data are given in Table 1. The 

data give the better determinations, both because of the larger cross- 
section and because of the reduced uncertainty from the PDFs when the 
better known u quark distribution is dominant. 

At HERA-II our ability to measure electroweak parameters will be greatly 
improved both due to increased statistics and due to the polarization of the 
beams, as detailed in the contribution of F. Metlica [17]. For example, with 
1/&- 1 and P(e~) = -0.7, the error achievable will be AM W ~ 0.055GeV, 
c.f. the PDG value AMw ~ 0.049GeV from time-like processes. 
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Experiment 


beam 


M w 


ZEUS 


e+ 


81.4 ± 2.7(stat) ± 2.0(sys) ± 3.0(PDF) 


HI 


e+ 


80.9 ± 3.3(stat) ± 1.7(sys) ± 3.7(PDF) 


ZEUS 


e~ 


80.3 ± 2.1(stat) ± 1.2(sys) ± l.O(PDF) 


HI 


e~ 


79.9 ± 2.2(stat) ± 0.9(sys) ± 2.1(PDF) 



Table 1. Values of M w extracted from ZEUS and HI CC data 



3. Recent NLO QCD fits 

A large amount of new data has become available during the past cou- 
ple of years, in particular the recent measurements of inclusive DIS cross 
sections in ep interactions by HI and ZEUS and the inclusive high-Er jet 
data by DO and CDF. The improved precision of the data led to a new 
generation of global NLO DGLAP QCD analyses, such as MRST01 [19] 
and CTEQ6 [20] presented at the workshop by R. Thorne [21] and W.K. 
Tung [22]. B. Reisert [23] and E. Tassi [24] presented QCD fits performed 
by HI [10] and ZEUS using their respective data supplemented by the data 
from fixed target experiments. Improved quality of the fits is achieved due 
to more precise data as well as due to a new level of sophistication in the 
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fitting technique including a full treatment of available experimental corre- 
lated systematic uncertainties. 

The MRST01, CTEQ6 and ZEUS PDFs are compared in Fig. 8, where 
the error band illustrated is that from the ZEUS standard (ZEUS-S) analy- 
sis. There is good agreement of all these PDFs within experimental uncer- 
tainties. 
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Fig. 8. Comparison of the ZEUS, MRST2001 and CTEQ6 PDFs. The error hand 
is that of the ZEUS standard (ZEUS-S) global fit. 



3.1. Uncertainty of parton distributions from QCD fits 

A new feature of the recent QCD analyses is a systematic and pragmatic 
treatment of the uncertainties of the parton distribution functions and their 
physical predictions. One of the problems of the uncertainty estimation for 
fits with many data sets is related to a certain degree of inconsistency of the 
latter. Usually the one sigma error of a parameter in a fit is determined by 
variation of x 2 by one unit from the minimum. Very often, however, this 
rule becomes unrealistic. This is demonstrated in Fig. 9, where the distances 
from x 2 -minima of individual data sets to the global minimum by far exceeds 
the range allowed by the Ax 2 = 1 criterion. It is not possible to simply 
drop "inconsistent" data sets, as then the partons in some regions would 
lose important constraints. On the other hand the level of "inconsistency" 
should be reflected in the uncertainties of the PDFs. This can by achieved 
by modification of the \ 2 tolerance criterion to Ax 2 = T 2 [20, 19, 25, 26] 
where T stands for a tolerance which should be estimated from the level 
of (in) consistency of the data sets used in each particular QCD fit. In the 
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CTEQ6 


Ax 2 = 


100 


a s (Mf)= 
a,(M%) = 


0.1165 ±0.0065(exp) 


ZEUS 


A * 2 eff 


= 50 


0.1166 ± 0.0049 (exp) 








±0.0018(modeZ) ± O.OOA(theory) 


MRST01 




= 20 


a s (M 2 z ) = 


0.1190 ± 0.002(exp) ± 0.003(theory) 


HI 


Ax 2 


= 1 


a s {Ml) = 


0.115 ±0.0017(exp) 
- o 0005 (model) ± 0.005 {theory) 



Table 2. Values of a s {M%) and its error from different NLO QCD fits with different 
error tolerances. 



CTEQ6 fit the tolerance was taken to be 10 (Ax 2 = 100), as shown by the 
horizontal lines in Fig. 9 (right). The choices for T 2 in the QCD fits are 
listed in Table 2 and range from 1 to 100. 
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Fig. 9. Partial x 2 for data sets in the MRST01 fit as function of a s (M§) (left). Dis- 
tance along a parameter combination (eigenvalue 4, CTEQ6) from the \ 2 -minimum 
of an individual data set to the global minimum (right). In the neighborhood of 
the global minimum a distance of 1 corresponds to ^X^iobai ~ 1- 



The values of the strong coupling constant a s (M z ) obtained in the fits 
are also given in Table 2. They are remarkably consistent. However, the 
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estimates of the experimental uncertainties on a s (M|) are different due to 
different judgements on the Ax 2 criterion. This is not a contradiction, all 
choices are legitimate and reflect different emphases in the fits. For example, 
HI [10] uses the canonical Ax 2 = 1 after careful consistency checks of the 
two data sets (HI and BCDMS fj,p) used in the fit. The relative uncertainty 
bands for the gluon distribution obtained in the CTEQ6 and ZEUS-S fits 
are shown in Fig. 10. They illustrate a reasonable consistency of judgement 
on the experimental errors of the gluon PDF. 
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Fig. 10. Relative gluon uncertainties as determined in the CTEQ6 fit (left) and the 
ZEUS-S QCD fit (right). 



Thus, there are reasonable approaches to how to treat experimental 
statistical and systematical errors, how to take into account model uncer- 
tainties such as charm or bottom masses, and how to account for incom- 
patibilities of data sets. It is not so easy to estimate theoretical uncer- 
tainties, this is explored further in [21] and in Sec. 5.4. The model un- 
certainty which comes from the choice of parametric forms for the PDFs 
at the input scale also merits further investigation. The latter issue was 
studied by CTEQ [20] and HI. B. Reisert [23] presented an investigation of 
the parameter space using general forms of MRST type parameterisations 
x ■ PDF = ax b (l — x) c (\ + d\fx + ex) for the PDFs of the gluon, quarks and 
anti-quarks at the input scale. Starting with the parameters o, b, c for each 
PDF, an additional parameter was considered only when its introduction 
improved x 2 by more than one unit. The uncertainty envelope in this study 
was defined as an overlap of the experimental and model error bands of fits 
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Fig. 11. Single differential cross sections da/dQ 2 for JVC and CC processes in the 
e ± p interactions (left). Grey error bands are envelopes of the QCD fit results 
corresponding to variations of parametric forms. The boxes in the left figure for 
JVC and CC at moderate Q 2 are zoomed and shown on the right. 



with x 2 < Xbest + y/ZHfof- Here Xbest corres P on ds to the fit with the best 
X 2 ■ This criterion is constructed by analogy with the statistical error of x 2 ■ 
It is somewhat arbitrary and used here only provisionally. The uncertainty 
bands in Fig. 11 show that an assumption on the parametric form for the 
PDFs at an input scale of Qo = 4 GeV 2 influences the predictions of the 
QCD fit even at moderate Q 2 in the case of highly integrated observables 
like da/dQ 2 . 

3.2. Results from QCD fits 

The HERA data are crucial for determining the low x sea and gluon 
shapes. The ZEUS-S sea and gluon distributions are compared in Fig 12 [24]. 
The gluon density is much larger than the sea density for Q 2 > 5 GeV 2 , 
but for lower Q 2 the sea density continues to rise at low x (consistent with 
the rise in F% down to low Q 2 mentioned in Sec 2.3), whereas the gluon 
density is suppressed. This could be a signal that the conventional DGLAP 
formulation of NLO-QCD is inadequate in this region. Fig. 12 also shows 
data at very low Q 2 compared to the ZEUS-S fit. Such a fit clearly fails 
for Q 2 < 0.65 GeV 2 , even when the conservative error bands on the fit are 
considered. 

For this generation of PDFs the data which most strongly determine 
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Fig. 12. Gluon and sea distributions for the ZEUS-S fit (left). F2 data at very low 
Q 2 compared to the ZEUS-S fit. (right). 



the valence distributions are the fixed target data. However new high Q 2 
HERA data can put some constraint on the valence distributions and this is 
important because these data do not suffer from the uncertainties associated 
with nuclear target corrections. 

Both the NC and CC high-Q 2 data are very well described by the global 
PDF fits. ZEUS has also made a special fit to ZEUS data alone (ZEUS- 
CD) including the new e~p 98/99 and the preliminary e + p 99/00 high Q 2 
data [24]. In this fit these additional data sets were used instead of the 
fixed-target data to constrain the valence distributions. Fig. 13 compares 
the valence distributions from the ZEUS-S global fit to those for the ZEUS-0 
fit. The level of precision of the ZEUS-0 fit is approaching that of the global 
fit and its precision is statistics limited rather than systematics limited, so 
that improvement can be expected with higher luminosity HERA-II data. 
The systematic precision of high- a; (x > 0.7) measurements at HERA-II can 
also be improved further as explored in the contribution of M. Helbich [27]. 

There are further advantages to using HERA data alone. In the ZEUS-0 
fit, the high-x (i-valence distribution is determined by the high-Q 2 e + p CC 
data. In contrast in the global fits it is strongly determined by the NMC 
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F2 /F$ data. It has been suggested that such measurements are subject to 
significant uncertainty from deuteron binding corrections [28]. The ZEUS-0 
extraction does not suffer this uncertainty. 

The ZEUS-0 fit was made using the same form of parton parametriza- 
tion as the ZEUS-S global fit. For the global fits, parametrization depen- 
dence is not severe since, for example, the valence shapes are strongly con- 
strained across all x by the CCFR xF% data, and the d — u distribution is 
constrained by D,p target data. However, if HERA data alone are used, 
then these constraints are lost and parametrization dependence can be sig- 
nificant, as discussed in Sec 3.1 [23]. At HERA-II the precision measurement 
of xF% across all x should considerably reduce this uncertainty, and elimi- 
nate the uncertainty from heavy target corrections which is unavoidable in 
the CCFR xF% measurement. 
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Fig. 13. Valence distributions for the ZEUS-S fit (left). Valence distributions for 
the ZEUS-0 fit (right). 
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Results relevant to the charm structure function and the photon struc- 
ture function were extensively covered in the Hadron Final States Working 
Group and will not be discussed again in detail here. Results from HERA 
on the charmed structure function were reviewed for the Structure Function 
Working Group by O. Behnke [29]. M. Prybycien [30] reviewed the status of 
photon structure function measurements at LEP.A. De Roeck [31] presented 
new OPAL data from two-photon processes and S. Maxfield [32] reviewed 
measurements of real and virtual photon structure. There was also some 
progress in developing new PDF sets for the photon, accounting for modern 
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data and correct heavy quark treatment, as presented by S. Albino [33] and 
P. Jankowski [34]. 

5. Theory 

At this workshop the main emphasis of the theoretical contributions 
was the different ways in which one can calculate structure functions and 
the regions of applicability of these different approaches. Essentially there 
were four alternative methods which were outlined, all of which have seen 
significant progress, or at least new results. These are: 

1. Lattice QCD. 

2. Saturation type effects/colour glass condensates. 

3. /c^-factorization. 

4. Collinear Factorization. 

There were also some other talks which do not fall into these general 
categories. D. Haidt presented a consistency check for DGLAP evolution 
[35], examining the partons extracted from the measured values of i*2 and 
di*2 jd In Q 2 and checking that these are consistent with the evolution equa- 
tions. A discrepancy is found at low x and Q 2 . T. Lastovicka demonstrated 
that a good fit to structure functions may be obtained using a parameteri- 
zation determined by assuming a self similar structure, i.e. using the fractal 
dimensions for the structure functions [36]. A. Kotikov presented a fit to 
high x data using cuts determined by the region of large systematic errors 
extracting, for example, a s {M 2 z ) = 0.1174 ± 0.0007(stat) ± 0.0019(sys) ± 
O.OOlO(norm) from a nonsinglet fit [37]. D. Timashkov also presented an 
analytic formula for structure functions for all x and Q 2 based on expres- 
sions in the limiting cases Q 2 — > 0, x — > 1 and x — > [38]. However, the 
summary is based on the above four alternative procedures. 

5.1. Lattice QCD 

There has been significant progress in this area, and we were given a sum- 
mary by S. Capitani [39]. It is not possible to compute structure functions 
directly on the lattice because the parton distributions are defined on the 
light cone, while lattice simulations are done in Euclidean space. However, 
one can use the Operator Product Expansion and calculate moments. The 
main effort has been in the calculation of 1st, 2nd and 3rd moments of nons- 
inglet distributions, both for unpolarized and polarized structure functions. 
A reason for only calculating nonsinglet quantities is due to the difficulty in 
computing disconnected diagrams (i.e., connected only by gluon lines) due 
to the expense in computer time. Nonsinglet quantities are insensitive to 
such diagrams. 
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One of the main improvements has been the first calculations without 
using the quenching approximation. This has shown that, for nonsinglet 
quantities at least, the quenched approximation is indeed very good. There 
have also been improvements in the perturbative renormalization factors 
required to translate the results on the lattice to a particular continuum 
renormalization scheme (e.g MS). In order to obtain the final results on the 
lattice it is ultimately necessary to perform chiral and continuum extrapo- 
lations, using a fit formula A + Bm 2 +ca 2 , due to the finite lattice spacing a 
and to the fact that one currently has a pion with mass ~ 500MeV. This 
appears to be well under control, but the results are disappointing - for the 
first moment of the u — d distribution they find 0.30 ± 0.03 where the stan- 
dard distributions give 0.23 ± 0.02. The results for polarized distributions 
are more in agreement with experiment. 

It is thought that this discrepancy is due to the finite size of the lattice 
missing the effects of the pion cloud. From chiral perturbation theory one 
obtain terms ~ m^lnmj. Introducing an additional term in the extrapola- 
tion formula of the form m 2 ln(m 2 / A 2 ) can solve this problem, but only for 
A a free parameter ~ 300 — 700MeV for various processes, destroying any 
predictive power. One needs lattices such that < 250MeV, which may 
be possible with computers within a couple of years. It would also be de- 
sirable to investigate the pion cloud effects by doing simulations with larger 
physical volumes. Finally, preliminary investigates of higher twist moments 
have been performed, giving results which are surprisingly small. 

5.2. Saturation type effects/colour glass condensates 

There was a lot of emphasis on the region of small x and low Q 2 where 
the gluon density is expected to saturate. Recently a great deal of interest 
has focused on the GBW saturation model [40]. In this one factorizes deep 
inelastic scattering into the fluctuation of the virtual photon into a dipolc 
pair and the dipole-proton cross-section. The former is calculable at LO 
and the latter is modelled in the form 

a{x,r) = <7 [1 -exp(-r 2 /4^(x))] 

where r is the dipole size and Rq(x) = (x/xo) x GeV~ 2 . This then saturates 
at large r and low x and predicts geometric scaling, i.e that the structure 
functions are functions of Q 2 Rq(x) alone. For A ~ 0.28 this model could 
be made to fit data well, and the procedure could then predict the total 
diffractive cross-section. We had a presentation from N. Timneanu on an 
extension of this type of approach to two photon physics [41]. Using a 
number of sensible extensions of the dipole-proton cross-section to model 
the dipole-dipole cross-section he demonstrated that a good fit to existing 
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data on photon-photon cross-sections is obtained for both real and virtual 
photons. However, this data is not very precise, and in particular, can not 
distinguish between the different extensions. 

However, the new more accurate HERA data on proton structure func- 
tions presented at this workshop is very precise and forces the simple satu- 
ration model to be modified, e.g., the effective power A in F2{x,Q 2 ) ~ x~ x 
is accurately measured to run with Q 2 (see Sec. 2.3 and Fig. 4) whereas it 
saturates at moderate Q 2 in the simplistic saturation model. In order to 
rectify this it is necessary to include DGLAP evolution in the model [42], 
replacing the above a(x,r) by 



<t[x, r 



1 — exp 



it r a s (fi )xg(x, /U 
3cr 



where fj, 2 is parameterized in terms of r 2 . This can improve the fit consid- 
erably, as shown in Fig. 14. However, it moves one away from the original 
principles somewhat, and geometric scaling is violated. There are two alter- 
native formulations in [42] (somewhere between the two is likely to be most 
realistic), both moving the saturation scale to lower x and Q 2 . Including 
charm quarks rather than only 3 light flavours in the original model (hardly 
optional since charm contributes over 30% of the structure function in some 
regions) also moved the saturation scale down an order of magnitude in x. 
If this is also true for the modified model it implies that for Q 2 = lGeV 2 
saturation sets in only for x < 0.00001 in structure functions. 




Fig. 14. The effective slope A as a function of Q 2 - the original saturation model 
(dashed line) and the improved model (dotted line). 



We also had some other talks on saturation in QCD, discussing the so- 
lutions to the nonlinear equations describing QCD at high parton density 
and adding support to the type of models considered above. We had a 
review of this subject by L. McLerran [43] (which will be discussed in the 
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summary of the diffractive session). We also had a presentation by A. Fre- 
und which outlined some specific points [44]. He told us that "Colour Glass 
Condensate" simply means that in QCD (colour) at high energies, fields 
slowly evolve relative to natural scales (glass) and the phase-space density 
saturates (condensate). Hence, at high energies, but at scales high enough 
that perturbative QCD applies, one can write a renormalization group equa- 
tion in x in terms of suitable variables [45] . This results in a Fokker-Plank 
equation which is nonlinear, but which one can use known techniques to 
solve. Implementing a single input condition one obtains a solution (for 
fixed as) of the general form F 2 (x,Q 2 ) = (x/x ) 2X f(x , ((^) A ^) 2 ) where 
A = 0.18, which fits the data well and is similar to geometric scaling. This 
has also been applied to nuclear structure functions with successful results, 
and where it is found that the scaling variable is Q 2 (x/xq) 2X A~ s , where 
5 = 0.1 rather than the naively expected 1/3. Hence, the first principle 
solutions in the colour-glass-condensate approach support the saturation 
models, implying a similar type of scaling. Corrections to this scaling have 
been estimated. However, there are further improvements to be made, e.g. 
a full treatment of running coupling, and it is intriguing that geometric 
scaling should work so well when the charm quark contribution is ignored, 
despite the fact that it contributes a great deal to the structure function, 
and should lead to violation of any geometric type scaling. 



5.3. kr factorization 

This should be applicable at high energies at scales where perturbation 
theory holds but high density effects are minimal. There have been a number 
of improvements in this field. One of these is in the Monte Carlos based on 
/cT-factorization. H. Jung showed that there has been significant progress 
in correcting previous shortcomings in the CASCADE Monte Carlo [46], 
both in the treatment of the scale in the running of the coupling and in 
the inclusion of the non-singular terms in the 0{as) gluon-gluon splitting 
function [47]. Explicitly, in the original version the splitting function had 
the form P = j^l + ^fA ns where A ns is the non-Sudakov form factor. This 
misses the non-singular terms as(— 2 + z(l — z)) in the LO splitting function. 
Since this is a negative contribution its omission leads to a bigger gluon at 
high and moderate x. This has lead to the modification 

P = a s {-^— z + Bz(l - z)j + a s + (1 - B)z(l - z)j A ns 

where B ~ 0.5, as well as changes to the form factors. Unfortunately, 
although these modifications, or something similar, are necessary, the de- 
creased positive contribution to the evolution at moderate x actually leads 
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to worse agreement with data for forward jets, and Tevatron b production, 
with the data showing an excess in both cases. There was also an alternative 
approach to fcy-factorization based Monte Carlos presented by G. Miu [48]. 
This is based on the Linked-Dipole-Chain model, and differs mainly in the 
manner in which partons are separated into initial and final state emissions. 
The resulting integrated gluon distribution obtained from fits to F2 agree 
well with standard distributions. While this need not be the case at small 
x, where /c^-factorization and collinear factorization may well differ, it must 
be the case at higher x, where a correctly modified Monte Carlo should not 
significantly alter the conventional results. 

A. Stasto talked on solutions to the LO BFKL equation with running 
coupling [49]. She argued that if one calculates the purely perturbative 
contribution to the high-energy gluon Green's function one obtains an ex- 
pansion in Poa'gY which is reasonably well-behaved as long as Poa s Y < 0.1 
(beyond this the series diverges), and explicit results are known for this 
series [50]. It was also demonstrated that the transition to the nonpertur- 
bative region is a sudden tunneling-like effect, rather than due to diffusion 
as is generally assumed [51]. The regime where the alternative methods of 
breakdown occur was compared and found to be similar for a>s < 0.1 but 
the perturbative expansion having a larger range of applicability for lower 
as (or in the formal limit of small (3$). However, f3oa s Y < 0.1 is not in 
practice a very wide range, and it is unclear if purely perturbative calcu- 
lations of high-energy scattering are really possible in a quantitative sense, 
though higher order corrections may help matters. 

S. Gieseke presented an update of the present status of the calculation 
of NLO impact factors in the BFKL framework [52]. This consists of two 
different contributions - the one-loop virtual corrections to the quark box 
diagrams and the contributions with an additional gluon in the intermediate 
state, as illustrated in Fig. 15. The calculation of the vertex diagrams in 
each of these two cases is now complete, and moreover, it has been proven 
that the infrared divergences due to the two separate contributions cancel 
each other in the appropriate manner. However, it still remains to perform 
the integrals over phase space to obtain the final result. Once this is done 
we will finally be in a position where NLO calculations of physical processes 
can be made with the fe^-factorization framework for the first time. 

5.4- Collinear factorization 

This is what we normally think of as standard perturbative QCD, which 
we use assuming high accuracy at large scales and uncertain x, usually at 
NLO. Within the framework of NLO QCD, S. Kretzer gave a presentation of 
various issues in heavy-flavour production [53]: a prescription for fully dif- 
ferential charm production at NLO for charm production in DIS [54], which 
will be important in determining the strange quark distribution with more 
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(a) (b) 

1 loop i tree tree i tree 




Fig. 15. Contributions to the 7* impact factor. 



accuracy at NuTeV; the NLO and mass corrections to the DIS contribution 
to u T N — > tX; and a modification of the ACOT prescription for inclusive 
charm production in DIS to include the appropriate threshold behaviour 
in coefficient functions at each order [55]. J. Bartels also gave a summary 
of the contribution of higher twist operators at small x [56]. Within the 
context of the double-leading-log approximation, and having to model the 
inputs, the contributions of the four sources due to gluon operators show 
a potentially large cancellation in F2, but imply a large negative twist four 
correction to F^. 

S. Moch gave a summary of NNLO calculations of splitting functions 
and coefficient functions [57]. These rely on calculating the Mellin mo- 
ments of the structure functions, which results in simplification since in- 
ternal propagators in diagrams which depend on the parton momentum p 
can be expanded in powers of (p • q/q 2 ) N where N indeed corresponds to 
the Mellin moment variable. For the diagrams with only one internal line 
dependent on p, known as basic building blocks, this reduces 4-point dia- 
grams to 2-point diagrams and the calculation is greatly simplified. Various 
techniques then also have to be used to relate more complicated diagrams 
to these building blocks [58]. A number of fixed moments of 3-loop splitting 
functions and coefficient functions have already been calculated [59]. The 
complete calculation of non-singlet quantities is nearly finished. The much 
more complicated singlet quantities will be a little longer. 

Finally, we had an update on both the MRST and CTEQ parton dis- 
tributions [19, 20]. W.K.Tung [22] concentrated on the treatment of uncer- 
tainties due to experimental errors (discussed in Sec. 3.1). R. Thorne [21] 
instead emphasized the need to understand theory errors for partons as well 
as the development of experimental errors. For example, MRST have used 
the approximate NNLO splitting functions in [62] to perform global fits 
and make predictions at NNLO [61]. This suggests that NNLO leads to a 
bigger correction even to the W cross-section at Tevatron than the experi- 
mental error within NLO. The theory errors associated with higher orders 
are probably much bigger for gluon dominated quantities. Additionally, de- 
tailed investigation of cuts on data [60] suggest that the fits improve if the 
lowest Q 2 and particularly lowest x data are cut out, and the predictions 
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for cross-sections with the new partons change. This suggests potentially- 
large corrections to NLO DGLAP at low Q 2 and low x. 

5.5. Theoretical Conclusions 

As outlined, there are various different approaches to calculating struc- 
ture functions, and there has been real progress in some of these calculations, 
e.g. NNLO in the usual expansion in as, in the colour glass model and in 
NLO corrections to BFKL impact factors. All the approaches are probably 
applicable in their own regime, and in some cases can be extrapolated with 
considerable success for surprising distances. However, this success may 
sometimes lead to unwarranted claims that one approach is actually partic- 
ularly appropriate. There needs to be more real understanding of precisely 
where the regimes are and how they can be combined in order to produce 
the best overall theory with the maximum predictive power. In our opinion 
the best theory, particularly when one considers the predictive power for a 
wide range of processes over a range of different experiments (HERA, Teva- 
tron, NuTev, LEP, LHC) is probably the collinear factorization theorem, 
but improved as much as possible by, for example, resummations at small 
and large x, higher twist corrections, etc. Clearly this will then require 
modification in the nonperturbative regime. This construction of the best 
complete, universally applicable theory is a difficult task, and help will be 
needed from even more precise and wide-ranging data. 

6. Conclusions 

New precise results on structure functions across a wide range of Q 2 , 
from 0.35 to 30000GeV 2 , have been presented this year. The data are now 
systematics rather than statistics limited such that QCD fits to extract 
parton distributions and a s have to consider correlated systematic errors. 
The precision of the data requires extension of the conventional formalism 
of NLO QCD, as embodied in the DGLAP formalism, and this challenge is 
being met as this formalism extended in various directions: to NNLO, to 
small x, to high density and to the non-perturbative regime. 
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